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SIGNAL  DETECTION  FOR  SAME-SHAPE  FAMILIES 


OF  NON-HOMOGENEOUS  POISSON  PROCESSES 

Roger  Lui*,  C.  B.  Bell,  Lloyd  Gavin**,  Joseph  Moser,  and  Edward  Pugh 


0.  Introduction  and  Summary 

Among  the  non-Gaussian  signal  detection  models,  the  models 
associated  with  non-homogeneous  Poisson  processes  (NHPP's)  are 
perhaps  the  most  applicable  and  most  tractable. 

This  paper  concerns  same-shape  families  of  NHPP's,  that  is, 
families  for  which  the  ratio  of  any  two  mean  functions  is  a  constant. 
One  such  family  is,  of  course,  the  family  of  homogeneous  Poisson 
processes  (HPP's).  Another  such  family  is  the  set  of  all  NHPP's 
with  mean  functions  of  the  form  u(t)  =  B  In  (l  +  3t),  where  B  >  0. 

In  general,  one  Is  concerned  with  n(PQ)»  which  contains  all  NHPP 
laws  with  mean  functions  of  the  form  p(t)  =  B  PQ(t),  where  B  >  0. 

For  each  such  family,  five  different  signal  detection  problems 
are  considered.  In  each  problem  the  pure  noise  (PN)  situation  is 
specified,  and  any  data  which  indicates  that  something  other  than 
the  specified  situation  has  occurred,  leads  to  the  conclusion  that 
some  signal  is  present,  i.e.,  one  has  a  noise-plus-signal  (N+S) 
situation. 

These  five  signal  detection  problems  are  as  follows. 

Problem  A  PN:  n(y  ) 

1  o 

(Goodness-of-fit  with  nuisance  parameter  or  class-fit) 
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Problem  B  PN:  fl(yQ)  and  u(.)  =  B0  yQ(  * ) 


( Goodness-of-f it ) 


Problem  C  PN:  =^f2  f£  n^o)l 


(2-sample  problem) 


Problem  D  PN:  ^  e  Q(y0)  and  • )  =  ot0  u2( ' ) 

(modified  2-sample  problem) 

Problem  E  <£.  e  ft(y0)  and  o^l  ~^-2  =  •••  =°C; 

(c-sample  problem) 


For  each  of  these  ( non-Gauss ian)  signal  detection  problems,  one 
considers  five  different  data  collection  schemes  or  sampling  plans. 

At  least  one  detection  rule  is  developed  for  every  one  of  these 
25  sampling-plan-problem  combinations. 

Section  1  introduces  the  notation,  distributions,  etc.  Section  2 
treats  the  sampling  plans  and  their  associated  likelihood  functions. 
Sections  3,  5,  6,  and  7  treat  detection  problems  A,  B,  C,  D  and  E, 

respectively.  Finally,  in  the  appendix,  numerical  illustrations  of 
each  technique  are  given. 
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Statistical  Preliminaries 


Most  of  the  distributions  and  concepts  employed  in  this  paper 
are  fairly  standard.  However,  it  is  worthwhile  to  give  the  more 
common  notation  and  abbreviations  one  will  encounter. 


1.1  Notation 

(1)  MSS  =  minimal  sufficient  statistic 

(2)  MLE  =  maximum  likelihood  estimate 

(3)  MVUE  =  minimum  variance  unbiased  estimate 

(h)  LHR  =  likelihood  ratio 

(5)  GOF  =  goodness-of-fit 

(6)  HPP  *  homogeneous  Poisson  Process 

(7)  NHPP  *  non-homo geneous  Poisson  Process 

(8)  PN  *  pure  noise 

(9)  N+S  ■  noise  plus  signal 

(10)  K-S  *  Kolmogorov-Smirnov 

(11)  i\-6-  Neyman-Barton 

(12)  0(uo)  *  the  family  of  laws  of  NHPP's  with  mean  functions 
of  the  form  »(•)  *  6  y0( • ) 


1.2  Methodology.  The  paper  concerns  NHPP's  and  inference  relative 
to  their  parameter  values.  The  situation  is  that  when  there  is  PN  (pure 
noise),  the  parameter  is  in  a  certain  range  of  values  and  when  there  is 
signal  present,  i.e.,  the  N+S  case,  the  parameter  is  in  some  other 
range  of  values. 
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This  being  the  case,  one  will  he  concerned  with  MSS's  (minimal 
sufficient  statistics)  for  the  parameters,  as  well  as  MLE's 
(maximum  likelihood  estimates)  and  MVUE's  (minimum  variance  unbiased 
estimate) . 

The  detection  procedures  will  be  based  on  reducing  the  detection 
problems  to  GOF  (goodness-of-fit),  2-sample,  and  c-sample  statistical 
problems.  The  GOF  problems  will  be  treated  with  K-S  (Kolmogorov- 
Smirnov)  and  (Neyman-Barton)  type  statistics.  (See  Section  1.3 
below)  For  the  2-sample  and  c-sample  problems,  the  detection  pro¬ 
cedures  will  be  LHR  (likelihood  ratio)  procedures,  or  conditional 
procedures . 

For  each  of  the  detection  problems,  the  number  of  possible 
procedures  is  quite  large.  In  this  paper,  only  one  or  two  procedures 
is  given  for  each  problem.  These  procedures  were  chosen  on  the  basis 
of  optimality  and  simplicity.  Some  of  the  references  in  the 
bibliography  contain  alternate  statistical  technicues  which  can  be 
adapted  to  signal  detection. 

1.3  Statistical  Distributions .  The  signal  detection  models 
considered  in  this  paper  involve  only  NHPP's.  However,  for  each  of 
the  detection  problems  five  different  sampling  plans  are  considered — 

Plan  A  (Type  I  Censoring); 

Plan  B  (Type  II  Censoring); 

Plan  C  ( Regular  Sampling ) ; 

Plan  D  (Equal-Distance  Sampling);  and 

Plan  E  (Same-Shape  Sampling). 


A  variety  of  standard  statistical  distributions  come  into  play. 


They  are  as  follows : 


(1)  Exp(A) 

(2)  P0( A) 

(3)  x? 


(M  Ge(p) 

(5)  U(O,0) 

(6)  r(n, A) 
(T)  B(m,p) 
(8)  F(m,k) 


(9)  N-B(k,p) 


(10)  €(•) 


=  Exponential  distribution  with  mean  a 

=  Poisson  distribution  with  mean  x 

*  Chi-square  distribution  with  "m"  degrees  of 
freedom 

=  Geometric  distribution  with  parameter  p 

=  Uniform  distribution  on  (0,0) 

=  Gamma  distribution  with  parameters  n  and  A 

=  Binomial  distribution  with  parameters  m  and  p 

=  Fisher's  F-distribution  with  m  and  k  degrees  of 
freedom,  respectively 

=  Negative  Binomial  distribution  with  parameters 
k  and  p 

*  the  one-point  distribution  with  all  probability 
at  0,  i.e.  e(u)  «*  1  for  u  >  D;  *  0  for  u  <  0. 


Besides  the  distributions  above,  one  will  be  concerned  with  several 
distributions  based  on  random  samples,  X^,  ....  Xm  i.i.d.  G(*),  with 
order  statistics  X(l),  ...,  X(m). 


(a)  G-O-S(m) 

(b)  U-O-S(m) 

(c)  EDF 


(d)  K-S(m) 


(e)K\-0(k) 


*  the  distribution  of  the  vector  [X(l),  ...»  x(m)] 

*  the  distribution  above  when  G  *  U(o,l) 

*  empirical  distribution  function, 

-1  m 

Fn(z)  =  m  2  e(z-Xj)  and  e(*)  is  as  in  (10)  above. 
This  is  a  mixture  of  one-point  distributions. 

=  the  distribution  of  the  Kolmogorov—Smirnov 

statistic  D(G,m)  *  sup  |  Fm(z)  -  G(z)|. 

z 

=  the  (exact)  distribution  of  the  Neyman- Barton 
statistic  (V-u)  £  ^  (V  -  u)',  where 


us  f—  -i-  — i— 1  v  =  fv^  Vkl 

~  2’  3*  *•••  k+1  ’  -  1  ’  ’  ***’  J 

with  V1*  *  k  ^  2  Xr,  and  2  =  (crs),  with 
*  rs  f (r+l)(s+l) (r+s+l) ]  1. 


2 

[This  distribution  is  asymptotically  x£.] 

For  each  of  these  distributions  one  writes,  e.g.,  X  ~  G,  to 
indicate  that  G  is  the  distribution  of  X;  and  [Y]_,  Ym]  ~  G-O-S(m), 

when  the  vector  has  the  distribution  in  (a)  above. 

One  now  introduces  the  five  sampling  plans  for  dealing  with  the 
NHPP's.  Several  of  these  plans  are  "natural."  For  example,  it  is 

often  necessary  (for  a  variety  of  reasons)  to  stop  sampling  at  some 

* 

fixed  time  T  ,  or  to  stop  after  the  occurrence  of  some  fixed  number 
of  events. 

Other  times,  it  is  desirable  to  sample  in  such  a  way  that  certain 
distributional  and  computational  difficulties  are  surmounted.  Such 
is  the  case  with  equal-distance  and  same-shape  sampling. 

All  of  these  plans,  their  associated  likelihood  functions,  MSS's 
and  MLE's  are  treated  in  the  next  section. 


Sampling  Plans,  Likelihood  Functions  MLE  and  MSS's 
Five  different  sampling  plans  will  be  considered. 


2.1  Sampling  Plan  A:  Type  I  Censoring.  The  stopping  rule 
here  is:  Stop  a  time  T*.  The  data  is  then, 

Z  =  [k,  Vlt  ...,  Wk] 

where  k  =  N(T*).  It  is  well  known  that 


Theorem  2.1.1  Conditionally,  given  N(T*)  =  k, 

(a)  W  =  [W]_,  ...,  Wk]  -  GPo-0-S(k), 

where  Gu(+)  =  ,  for  0  <  t  <  T*,  and 

u(T*) 

^wy0(wl)  Uo(Wk) 

(h)  [  *»>  •••»  ,  # .  J  “  UOS(k). 

U0(T*)  u0(t*) 

(c)  Further,  the  likelihood  function  is 

L(z)  =  L(w  |  N(T*)  =  k)  P{  N(T#)  =  k>  =  gk  [k  v'0  (wj)]  exp  {  -  g  ^(T*)} 

1 


From  the  likelihood  function,  one  derives  the  following. 

Theorem  2.1.2  (a)  The  MSS  for  g  is  N(T*)  -  PQ( guo(T*) ) . 

(b)  The  MLE  of  g  is  g  =  and 

(c)  E(g)  =  g  and  V(  g)  =  — -  . 

Uo(T  ) 

The  signal  detection  methodology  for  this  sampling  plan  will  employ 
the  previous  two  theorems. 


2.2 


iling  Plan  B: 


»S3= 


II  Censorim 


-  8  - 


£ 

Theorem  2.2.1.  (a)  Conditionally,  given  Wn  =  w  , 

W*  =  [WX,  %_i]  ~  G^-O-SCN-l),  and 

r^Wl)  Wo(Wr,_n  )  . . .  ,, 


(b)  [■ 


'T7^’  “oW 


-]  ~  U-O-S(N-l). 


(c)  Further,  the  likelihood  function  is 
L  (w)  =  L  (vi,  ...,  w]j_i  |  Wjf  =  w*)  fWN  (w*) 

■  sN  fn  u'  (wj)]  exP  {-eu0(v*)}  . 


Employing  the  likelihood  function,  one  obtains 


Theorem  2.2.2  (a)  The  MSS  of  6  is  u0(Wn)  **  r(N,B). 
Further,  2Bv0(WN)  -  X2 


(h)  The  MLE  of  B  is  B 


Mo(%) 


with  E(b)  = 


and  V(b) 


(N-l)2(N-2) 


(c)  The  UMVUE  of  B  is  8*  =  — SlL  =  (Ili)  g,  with 

VoW 

E(B*)  =  B  and  V(8*)  = 


The  signal  detection  procedures  for  this  sampling  plan  will 
primarily  be  based  on  the  two  preceding  theorems. 


Lng  Plan  C: 


Lng.  Here,  one  observes 


the  process  at  times  A,  2A,  ...,  kA,  and  the  data  is  Y  =  [Y]_,  ...»  Yk] 

where  Yi  *  N(a),  Y2  =  N(2A)  -  N(a),  ...  Yk  *  N(kA )  -  N([k-l]A). 


The  basic  distribution  results  are  as  follows: 


Theorem  2.3.1  (a)  Y^,  ...,  Yk  are  independent. 

(b)  Yj  -  P0(bj6),  where  bj  =  u0(jA)  -  P0([J“l]A)»  for  j  =  1,  2,  . k. 


(c)  The  likelihood  function  is 

L(y)  =  B1^1^  [  *  bjy^  ]  [  *  yj!]  exp  {-  8y0(kA )} 

1  1 

One  uses  the  likelihood  function  to  establish 

Theorem  2.3.2  (a)  The  MSS  for  6  is  N(kA)  -  P0(Bu0(kA)) 

(b)  The  MLE  of  8  is  8  =  ^T£a) 

(c)  E(B)  =  6  and  V(8)  = 

As  is  usual,  these  last  two  theorems  form  the  basis  for  the 
signal  detection  procedures  for  regular  sampling. 

2.h  Sampling  Plan  D:  Equal-Distance  Sampling.  In  this  sampling 
plan  one  observes  N(ti),  N(t2),  ...,  N(tk)  where  pD(tr)  =  r  u0(ti) 

The  data  here  is  Y  =  [Y]_,  ...,  YkJ,  where  Yr  =  N(tr)  -  N(tr_i) 
for  r  =  1,  2,  ...,  k  and  tQ  =  0. 

The  pertinent  distribution  result  is 

Theorem  2.4.1  (a)  Y^,  ...,  Yk  are  i.i.d.  P^Buq^i)) 

(b)  The  likelihood  function  is 

L(y)  =  [8u0(ti) ]N(tfc)  [  n  yj!]  exp  {-  kBu0(ti)}. 

Immediate  consequences  are  as  follows. 


Theorem  2.^.2  (r)  The  MSS  for  6  is  N(tk)  -  PQ(  Buo^k)  )■ 

,  .  -  N(tk)  N(tk) 

(b)  The  MLE  of  6  is  e  =  ^ 


s  , 

\\ 


(c)  E(b)  =  8  and  V(b)  =  — j \ 


These  last  two  theorems  will  be  utilized  in  constructing  the 
detection  procedures  for  the  equal-distance  sampling  plan. 

The  final  sampling  plan  is  the  same-shape  sampling  plan  adapted 
from  Basawa  and  Rao  (1980),  who  were  concerned  with  homogeneous 
Poisson  Processes  (HPP's). 

2.5  Sampling  Plan  E:  Same-Shape  Sampling.  Let  (N(t)}  ,  the 
NHPP  of  interest,  have  mean  function  By0( • )  and  waiting  times  {Wr}, 
and  let  {N*(t ) }  be  an  NHPP  independent  of  (N(t)}  and  having  mean 
function  u0( • )  and  waiting  times  {W*}  . 

One  observes  N*  =  [N(W*),  ...»  N(W*)]  and  computes  Y  =  [Y]_,  ... 

Y1  =  N(W*),  Y2  =  N{W*)  -  N(W*),  . . . ,  Yfc  =  N(W*)  -  NfwJ^)  . 

It  can  be  readily  proved  that 

Theorem  2.5.1  (a)  Y^,  ...,  Yfc  are  i.i.d.  Ge(p)  with  p  =  [l  +  6] 

(b)  |  Y  =  N(W*)  -  N-B(k,p). 

1  K 

(c)  The  likelihood  function 

Ky)  =  pk  /'“k1  *  SH(Mk>  [1  .  6]-k-"(WJ) 

It  follows  that 

Theorem  2.5.2  (a)  The  MSS  for  6  is  N(W*). 

„  jf  ( ) 

(b)  The  MLE  of  g  is  8  =  - —  =  Y. 

k 

.  (c)  E(S)  =  B  and  V(B)  — 

One  is  now  in  a  position  to  develop  the  various  signal  detection 


procedures . 


The  detection  rules  will  be  designates-  first,  by  problem;  second 
by  sampling  plan;  and  third,  by  the  number .  For  example,  "Detection 
Rule  B.II.l"  refers  to  the  first  detection  rule  for  Detection 
Problem  B  with  Sampling  Plan  II;  and  "Detection  Rule  C.III.2"  refers 
to  the  second  detection  rule  for  Problem  C  and  Sampling  Plan  III. 


3.  Signal  Detection  for  Problem  A:  Class-fit,  or  Goodness-of-fit 


with  Nuisance  Parameter. 

The  situation  here  is  that,  if  there  is  pure  noise,  the  NHPP  law 
is  in  while  if  signal  if  present,  what  is  received  is  governed 

by  a  law  not  in  n(vQ).  One  has 

PN  :  £  e  fl(Hp)  vs  N+S_  ■  • 

The  decision  procedures,  of  course,  depend  on  the  sampling  plain. 


3.1  Sampling  Plan  A  (Type  I  Censoring)  Here  one 
where  N(t*)  =  k.  and  computes  V  = [V^,  ...,  V^]  where 


stops  at  time  T*, 


Vr  = 


Wp(Wr) 
Uo(T*)  ' 


Then,  as  in  Theorem  2.1.1,  V  U-O-S(k). 

-  W 

Hence,  the  class-fit  problem  is  now  a  GOF  problem. 

One  now  forms,  F*(z)  =  —  £  e(z  -  V  )  ;  Vs  *  -r-  '  £  V?,  and 

k  k  i  r  K  ]_  <J 

v  =  [v1,  v2,  ...,  vPj  ,  and  bases  the  decision  rules  on  these  entities. 


Detection  Rule  A. 1,1  Decide  N  +  S  iff  D^  =  sup  |  f£  (z)  -  z 

o<z<l 

where  d^k  is  the  appropriate  percentile  from  a  K-S(k)  distribution. 

This  K-S  detection  rule  is  considered  an  "omnibus"  rule  in  the  sense 
that  it  is  consistent  against  "all"  alternatives. 

Adaptation  of  the  work  of  Neyman  (1937)  and  Barton  (1953)  yield 
procedures  recommended  for  certain  "smooth"  N  +  S  situations. 


Detection  Rule  A. 1.2  Decide  N  +  S  iff 
Bp  =  k  (V  -  v)  I"1  (V  -  y)'  >  Xp(l-o),  where  V  =  [j,  -j,  ...,  and 

^  ~  ^°rs^’  ars  (r+1)  (s+1)  (r+s+1) 


The  detection  procedures  are  different  for  different  sampling 
plains.  However,  for  this  detection  problem,  they  are  similar  for 
the  first  two  sampling  plans. 

3.2  Sampling  Plan  II  (Type  II  Censoring)  The  data  here  is 
W  *  [W15  ....  Wjj],  and  under  PN  [V*,  ...,  vj^]  ~  U-O-S  (N-l), 

*  po(Wr) 
where  Vr  =  — rrr  . 

V  V 

The  detection  rules  here  are 


Detection  Rule  A.II.l  Same  as  Detection  Rule  3.1.1  where 
k  is  replaced  by  N-l  and  V's  are  replaced  by  V  's. 


Detection  Rule  A. II. 2  Same  as  Detection  Rule  3.1.2,  where 
V's  are  replaced  by  V*'q. 


3.3  Sampling  Plan  III  (Regular  Sampling)  For  this  sampling 
plan,  Y^,  ...,  Yfc  are  independent  with  Yj  ~  PQ(bjB)  under  PN  (See 
Section  2.3). 

A  useful  straight forward  detection  rule  is 


Detection  Rule  A.III.l  Decide  N  +  S  iff 

k  [y  -=^“]2 
T  =  £  LyJ  b  •  = 

J=1  b.  __ 


[Wo(kA)  N(kA)]-1  £  bj1  [yj  v0(kA)  -  bj  N ( kA ) ] 2  >  x^-i  (a) ■ 


For  the  next  sampling  plan,  the  increments  yj  are  i.i.d.,  and,  hence 
the  detection  procedure  is  somewhat  simpler. 


3.1+  Sampling  Plan  IV  (Equal- Distance  Sampling)  One  recalls 


(Section  2.1+)  that  Y^  . ..,  Yk  are  i.i.d.  P0(8y0(ti)),  where 
Yr  =  N(tr)  -  N(tr_i),  t0  =  0,  and  yQ(tr)  =  ryQ(ti).  A 
decision  rule  of  interest  is  then 


Detection  Rule  A.IV.1  Decide  N  +  S  iff 

k 

X2  ( 1-  a)  <  T*  =  (Y)_1  E  (Y  -  Y)2 

k-1  j=l  J 


ij7T-T  E  t»(tr)  -  K(tr-i) 
' uk'  r=l 


W(tk)  2 
k 


Each  of  these  rules  is,  of  course,  illustrated  by  a  numerical 
example  in  the  appendix. 

The  final  detection  rule  for  this  section  is  based  on  same-shape 
sampling. 


3.5  Sampling  Plan  V  (Same-Shape  Sampling)  For  this  sampling 
plan,  under  PH  the  Y's  are  i.i.d.  Ge(p)  where  p  *  [1+0]  1,  and 
Yr  -  R(W*)  -  N(W*_!)  with  W*  =  0.  (See  Section  2.5) 

Following  usual  procedures  one  develops 


Detection  Rule  A.V.l 

[nr  -  kp(Cr)]2  > 
r=1  kp(Cr) 


Decide  N  +  S  iff 


( 1-e ) 


Here,  {A^,  ...,  Ag}  are  integers  satisfying  Aq  =  0  <  A^  <  A,, 

<  ...  <  A  <  «.  Then,  one  chooses  C  =  {1+A  - ,  2  +  A  A 

s  r  r-1  r-1  : 

for  1  <.  r  <  s,  and  C1*{A,2  +  A,...}. 


Let  n  =  the  number  of  Y's  in  C  ,  and 
r  r 


v  1+A  1 

p  (c  )  =  (-2L)  r_1  -  (-i=) 

r  l+Y  1+Y 


Y  1+Ar- 

1  '  r  .  for  1  <  r  <  s 


„  ^7  1+A  — 

and  P(C  )  =  1  +  (— )  r  -  (— )  . 

r+1  '  l+Y  1+Y 


{C^,  . ..,  constitute  a  partition  of  the  non-negative 

integers,  and  there  are  many  such  partitions.  The  partition  which 
will  be  used  in  the  numerical  examples  of  the  appendix  is  as  follows 

=  {0,5,10,15,...}  ;  =  {1,6,11,16,...}  ;  =  {2,7,12,17,...}  ; 

C£  =  {3,8,13,18,...}  and  =  {fc,9,lfc,19, . . . }  . 

>•—1  C  “1 

Here,  p(Cr)  =  p(q)  [l  -  ( q.) ^ ]  for  r  =  1,  ...,  5,  where 
p  =  [l+Y]  and  q  =  1  -  p. 


One  turns  now  to  the  second  detection  problem. 


h.  Signal  Detection  for  Problem  B:  Goodness-of-Fit 

The  PN  case  here  entails  the  stochastic  process  law  to  be  completely 
known,  i.e.  to  be  a  NHPP  law  with  mean  function  v(*)  satisfying 
p(t)  =  given  that  e  G(u0). 

PN  :  u( • )  =  Bopq( *  )  vs .  N+S  :  u( • )  ^  ^o^o^  * \ 

The  problem  here  is  to  decide  whether  or  not  8  -  BQ. 

l+.l  Sampling  Plan  I  (Type  I  Censoring)  Under  PN  (Section  2.1)  the 
MSS  for  8  is  N(T*)“  PQ( 80Po(T#^ •  ^  detection  rule,  then  is 

Detection  Rule  B.I.l  Decide  N  +  S  iff  N(T*)  _<  a1  or  N(T*)  >  a2, 

where  E1  - — p-  +  z  ~rT'~  =  °  with  X  =  Bou°^T  ^ 

r=0  •  r‘  rs*a2 

In  practice,  one  would  probably  wish  to  choose  a1  and  a^  to  have 
approximately  equal  tail  probabilities  under  PN. 

Contrary  to  what  was  true  for  Problem  A,  the  detection  rule  for 
Sampling  Plan  II  here  is  somewhat  different  from  that  for  Plan  I. 

k.2  Sampling  Plan  II  (Type  II  Censoring)  One  vises  the  fact  that 

o 

under  PN,  2B0yo(%)  *  XgN>  to  construct  a  reasonable  detection  rule. 

(See  Section  2.2) 

Detection  Rule  B.II.l  Decide  N  +  S  iff  280Vo(%)  <  b^  or  >  bg 

2 

where  b^  and  b^  are  appropriate  percentiles  of  the  X2JJ  “  distribution. 
For  the  third  sampling  plan,  one  makes  use  of  an  approximate 


chi-square  distribution. 


is  Y 


=  (Y  ,  ....  Y,  ),  where  Y  =  N(rA)  -  N([r-l]A),  with  N(0)  =  0. 

X  iC  I* 

Further,  one  recalls  (Section  2.3)  that  under  PN  Y, ,  . . . ,  Y,  are 

1’  ’  k 

independent  with  Y  -  P  (Bbj.),  where  bj.  =  u0(rA)  -  y0([r-l]A). 

One  useful  detection  rule  is 

Detection  Rule  B.III.l  Decide  N  +  S  iff 

T  =  E  (B0  br)"1  Hr  ~  60  ^]2  >  X?  (l  -  «)•. 
r=l  K 

Another  such  rule  (which  employs  the  MSS  for  B)  is 

Detection  Rule  B.III.2  Decide  N  +  S  iff  E  Y  =  N(kA)  <  an  or  >  a0 
-  !  r  .  "  1  ~  2 

(as  in  Detection  Rule  B.I.l). 

The  .situation  is  somewhat  different  for  equal-distance  sampling. 

k.U  Sampling  Plan  IV  (Equal-Distance  Sampling)  From  Section  4.U, 
one  recalls  that  yQ(tr)  =  rt^;  Yr  =  R(tr)  -  N(tr_i)  with  N(0)  =  0; 
and  that  Y^  . ..,  Yfc  are  i.i.d.  PQ(  Bo^o^) )  under  PN,  and  that 
N(t^)  is  the  MSS  for  B. 

Under  such  circumstances  one  constructs 

k 

Detection  Rule  B.IV.l  Decide  N  +  S  iff  E  Yr  =  N(tk)  £  a1  or  >  &2 
(as  in  Detection  Rule  B.III.2). 

^.5  Sampling  Plan  V  (Same-Shape  Sampling)  The  M-S-S  for  B  (Section  2.5) 

k  *  -1 
E  Y  =  N(Wk),  which  under  PN,  -  N-B(k,p)  for  p  =  [l  +  B0]  • 


Under  these  circumstances  one  employs  the  detection  rule  below. 


Detection  Rule  B.V.l  Decide  N  +  S  iff  I  Y  =  N(W*)  <  b,  or  >  b„. 

where  b1  and  b^  are  appropriate  percentiles  of  the  Negative  Binomial 
distribution. 

This  final  detection  rule  as  well  as  the  other  rules  of  this  section 
are  illustrated  in  the  appendix. 

For  the  next  detection  problem,  the  rules  are  quite  different  in 
the  sense  that  most  of  them  sure  conditional  rules. 


For  this  problem  it  is  assumed  that  data  from  two  processes  with 
laws  in  D(u0)  Eire  being  received.  The  PN  situation  entails  the  equality 
of  the  two  laws,  which  is  equivalent  to  the  equality  of  the  mean  functions 
8-jUo  and  82p0- 

PN  =^2  (i.e.  61  =  B2)  vs  N  +  s  :<^i  ^2 

5.1  Sampling  Plan  I  (Type  I  Censoring)  Let  {N(t)>  and  (M(t)}  be 
the  two  independent  NHPP  processes  and  {W^}  and  {W^ }  be  their  respective 

waiting  times.  Observation  of  both  processes  is  in  the  time  interval 
[0,T*] .  The  data  is  then  (following  Section  2.5) 

Z  =  [n;  W1,  . ..,  Wn;  m;  W*,  ....  ¥*] 

where  n  =  N(T*),  m  =  M(T#). 

N(T*)  and  M(T*)  are  the  M-S-S's  for  8^  and  82,  respectively.  They 
are  independent  with  N(T*)  *  )  an<*  M(T#)  *  P0(  B2jj0(t*))  • 

Hence, 

Lemma  5.1.1  Under  PN,  and  conditionally,  given  N(T#)  +  M(T*)  =  N, 
N(T*)  -  B(n,Js). 

The  "natural"  detection  rule  is  then 

Detection  Rule  C.I.l  Decide  N  +  S  iff  N(T*)  £  c^  or  >  c2,  where 
and  c2  are  the  appropriate  percentiles  of  a  B(N,Js)  distribution. 

This  detection  rule  is,  as  usual,  illustrated  by  a  numerical  example 


in  the  appendix. 


For  the  second  sampling  plan,  the  detection  rule  is  based  on  an 
F-distribution. 

5.2  Sampling  Plan  II  (Type  II  Censoring)  For  this  plan  the  data  is 

Z  =  [W, ,  W  ;  W*,  ...»  W*],  i.e.  one  observes  the  first  n  waiting 

-  x  n  l  m 

times  of  (N(t)}  and  the  first  m  waiting  times  of  (M(t)}.  Then,  one  has 

mu  (  W  ) 

Lemma  5.2.1  Under  PN,  —  -  F(2n.2m). 

ny0(w*) 

The  "natural"  detection  rule  here  is  then 

my0(Wn) 

Detection  Rules  C.II.l  Decide  N  +  S  iff - sr  <  f„  or  >  f^.  where 

-  nu(W*)  1  2 

o  m 

f^  and  f^  are  appropriate  percentiles  of  an  F(2n,2m)  -  distribution. 

For  the  next  two  sampling  plans , • the  detection  rules  are  based  on 
conditional  binomial  distributions  as  in  Section  5.1  above. 

5.3  Sampling  Plan  III  (Regular  Sampling)  The  M-S-S's  for  the 
6^  and  B2  8X6  N(nA^)  and  M(mA2),  respectively  (Section  2.3). 

Further,  it  is  immediate  that 

Lemma  5.3.1  Conditionally,  given  N(nA1)  +  M(mAg)  =  N, 

u0(nA1 ) 

N(nA  )  -  B(N,p),  where  p  - - 

1  F  U0(nA1)  +  uQ(mA2) 

[Note:  If  A^  =  A2  and  m  =  n,  then  p  =  h.] 

The  detection  rule  is  then 

Detection  Rule  C.III.l  Decide  N  +  S  iff  N^A^  <  a^  or  >  a2 
(where  a^  and  a2  are  appropriate  percentiles  of  the  B(N,p)  -  distribution) 


-  21  - 


An  analysis  rule  is  derived  for  equal-distance  sampling. 

5.1*  Sampling  Plan  IV  (Equal-Distance  Sampling)  The  data  here  is 


Z  =  [N(t-),  ....  N(t  )  M(t.),  ...,  M(t  )J  and  the  M-S-S  for  B.  and  B 

„  l  ni  m  -L 

are  H(t  )  and  M(t  ),  respectively, 
n  m 

The  detection  mile  is,  then, 


Detection  Rule  C.IV.l  Decide  N  +  S  iff  N(t_)  <  b..  or  >  tu  (where 


h.  and  b0  are  appropriate  percentiles  of  the  B(n,p)  -  distribution), 

u  (t  ) 

with  N  =  N(t  )  +  M(t  )  and  p  =  u  U  7  • 


u  (t  )  +  u  <t  ) 
on  o  m 


The  final  detection  rule  for  this  signal  detection  problem  is  also 
based  on  a  conditional  distribution.  However,  this  conditional  distribution 
is  not  binomial. 


V  (Same 


Lng)  The  two  independent 


stochastic  processes  of  interest  are  (N(t)}  and  (M(t)}  with  mean  functions 
g^u0( • )  and  BgUof'K  respectively.  One  introduces  two  independent 
NHPP's  (H*(t)}  and  {M*(t)J,  which  are  independent  of  the  original  processes, 
and  which  have  waiting  times  {W*>  and  {V*>,  respectively. 

The  data  here  is  then 

Z  =  [N(WJ),  ...,  N(W*);  M(vJ),  ...,  M(V*)] 


Prom  Section  2.5  and  some  straightforward  derivations,  one  concludes 

Theorem  C.5.1  (i)  N(W#)  and  M(V*)  are  the  M-S-S's  for  B-.  and  8„, 

n  m  c 

respectively,  (ii)  Under  PN,  P{  N(w#)  *  s  I  N(W*)  +  M(V*)  =  k  } 

n  n  m 

/  n  +  s-1  \  /  zn  +  k-s-1  \r/n  +  k“*l\i  0  /  i  .  \  ,  ,T 

*  (  _  )  (  v__  )  l(  .  )J  =  P(s  |  n,  m,  k)  where  N  =  n  +  m. 


-  22  - 

On  the  basis  of  this  result,  the  detection  rule  becomes 

Detection  Rule  C.V.l  Decide  N  +  S  iff  N(W*)  <  ^  or  >  c2,  where 
C1 

Z  P(s  |  n,  m,  k)  +  Z  P  (s  |  n,  m,  k)  = 
s=0  s=c2 


a. 


6.  Problem  D  (Constant  Intensity  Ratio) 

This  detection  problem  is  closely  related  to  the  previous  problem. 
Essentially,  for  Problem  C,  PN  entailed  the  equality  of  the  intensities 
of  two  independent  NHPP's.  Here,  PN  entails  that  the  intensity  of  the 
first  process  be  a  fixed  constant  multiple  of  the  intensity  of  the  second 
process. 

Under  such  circumstances,  the  data  vectors  and  the  M-S-S's  for  the 
various  sampling  plans  are  the  same  as  in  Section  5;  and  the  detection 
rules  are  closely  related  to  those  of  Section  5.  For  that  reason,  these 
rules  will  be  listed  with  a  minimum  of  discourse. 

PN  :  61  =  y062  ^  N  +  S  :  B1  i-  y0&2 


6.1  Sampling  Plan  I  (Type  I  Censoring) 

Detection  Rule  D.I.l  (See  Rule  C.I.l)  Decide  N  +  S  iff 

N(T*)  _<  or  _>  a g,  where  a^  and  a^  are  appropriate  percentiles  of 

Y 

the  B(N,p)  -  distribution,  with  N  =  N(T*)  +  M(T*)  and  p  =  -  -+°— ■ 

o 

6.2  Sampling  Plan  II  (Type  II  Censoring) 


Detection  Rule  D.II.l  (See  Rule  C.II.l)  Decide  N  +  S  iff 

m YoWo^V  ) 

- <  f,  or  >  f_,  where  the  f's  are  the  appropriate  percentiles 

“  «.(<> 


of  the  F(2n,2m)  -  distribution. 


6.3  Sampling  Plan  III  (Regular  Sampling) 


Detection  Rule  D.III.l  (See  Rule  C.III.l)  Decide  N  +  S  iff 
N(nA1)  <.  or  >  a2,  where  a1  and  a2  are  appropriate  percentiles  of 


the  B(N,p)  -  distribution,  with  N  =  n  +  w  and  p  = 


Youo(nA1) 


YoUo(nAl)  +  ^0^2 


6.U  Sampling  Plan  IV  (See  Rule  C.IV.l)  (Equal-Distance  Sampling) 


Detection  Rule  D.IV.l  Decide  N  +  S  iff  N(t  )  <  b,  or  >  b 


where  b^  and  b2  are  appropriate  percentiles  of  a  B(N,p)  -  distribution 

Y.u  (t  ) 

with  N  =  N(t  )  +  M(t  )  and  p  =  - 7— — r— - rr — ?  . 

n  m  Y„yJtJ  +  U  (t  ) 

o  o  n  o  m 


6.5  Sampling  Plan  V  (See  Rule  C.V.l)  (Same-Shape  Sampling)  The 


detection  rule  for  this  sampling  plan  is  quite  different  from  the  other 

rules  encountered  up  to  this  point.  The  crucial  point  here  is  that  the 

M-S-S's  N(W*)  and  N(V*)  are  distributed  N-B  (n,  _  J— ■)  and  N-B(m,  ■9~\ 

m  X  +  p.^  X  +  Yqo 

respectively  under  PN.  Hence,  their  sum,  when  yq  f  1»  is  not  distributed 
as  a  Negative  Binomial  variable.  One,  hence,  resorts  to  a  likelihood 
ratio  procedure. 

For  the  special  case  m  =  n  =  k,  one  has 


Theorem  6.5.1  (i)  The  likelihood  function  is 

/  \  /  ik  N(w£)  M(v£) 

L(yn,  ....  ylk;  y21,  •••,  72k>  =  (p1p2)  *1  q2 


where  p. 


and  q_  =  1  -  p  .  for  r  =  1,2. 


(ii)  The  unconditional  maximum  of  L(y)  is  achieved  for 


n(w£) 


M(VJ)  _! 


Px  =  f  1  +  — j —  ]  and  p2  =  [  1  +  — - — ] 


(iii)  Under  PN,  p,  =  - - - rs - r:  and  L(y)  is  maximized 

1  YQ  +  P2  (1“Y0)  - 


A  A  M 1  A 

for  p2  =  [l  +  x]  and  p^  = 


Y0  +  P2  (1_Y0) 


where  x  =  [U  y  ]_x  ({ [ 1  +  y  -  k"1  N(wJ)]2  +  8  y  [N(W*)  +  M(V*)]} 


1  f  /  T  T^  > 


-  [1  +  yq  -  k_1  N(W*) ]) 

.  .  H(WJ)  -k-N(w£)  M(V*)k-M(v*) 

(iv)  L(y,  p, ,  p 0)  =  [1  + - J  [1  + - J 

-  x  k  k 


k  k 


(v)  L(y,  Pj,  p2)  =  {yq  p2  (1-yo)  [(l-v0)  p2  +  Yq]  ) 

!p2  d-v0)]M<V^  [<1-T  ,  ;2  *  r  fk 


-a  »<«;> 


From  standard  statistical  derivations,  it  follows  that 

_  (y»  p-,,  p2)  ? 

Theorem  6.5-2  under  PN,  S  =  -  2  In  — = — : - : —  z  yT 


(y»  Pj_t  p2) 


Therefore,  the  detection  rule  is 


Detection  Rule  D.V.l  Decide  N  +  S  iff  S*  <  y2  (a). 


Signal  Detection  for  Problem  E:  c -Sample 


For  this  problem,  one  has  data  from  c  independent  NHPP's  : 

(N^(t)},  { N^( t ) } ,  (Nc(t)}.  The  laws  of  these  processes, 

...,^c  are  all  in  n(y0). 

The  detection  problem  is 

PH  :  ^  -  ...  vs  N  +  S  :  ' s  are  not  all  equal 

(i.e.  8X  =  $2  =  . . .  =  Bc) 

For  this  problem,  the  detection  rules  are  much  more  tractable  if  the 
sample  "sizes"  are  "equal,"  in  the  appropriate  sense,  and  this  latter 
condition  will  be  assumed. 

The  data  for  the  various  sampling  plans  is  as  given: 

Data  for  Sampling  Plan  I:  W  =  [W.^,  W,,,  . ,  WQ],  where  the  waiting 

times  are  =  (Wrl,  ...,  W^)  and  Nr(T*)  =  k  for  r  =  1,  2,  ...,  c, 

£ 

i.e.  each  process  is  observed  until  time  T  . 

Data  for  Sampling  Plan  II:  W  =  [W„ . W  ],  where  W  =  [W  W 

„  „x  „c  „r  rx  rN 

for  r  =  1,  2,  ...,  c;  I.e.  each  process  is  observed  until  the  Nth  waiting  time. 


Data  for  Sampling  Plan 


III:  N  =  [N1,  ...,  Nj  or  Y  =  [Yr  ...,  Yc] 


where  ^  =  [y*),  H.(2A),  ...,  N.(kA)]  and  Yr  =  fYrl,  ....  Yrk], 
with  Yrg  =  N^sd)  -  N  ( [s-l]d)  and  N  (0)  =  0  for  r  =  1,  2,  ...,  c 


1,  2,  •  •  • , 


/\ 


Data  for  Sampling  Plan  IV:  H  =  [N^,  or  Y  =  [Y^,  Y^j 

where  Nr  =  [N^t^,  . ..,  Nr(tfc)]  and  Yr  =  [Yrl,  ....  Yrk],  with 

Y  =N(t)-N(t  t)  and  N  ( 0 )  =  0  for  r  =  I ,  2 ,  . . . ,  c  and 

rs  r  s  r  s-l  r 

s  =  1,  2,  . . . ,  k. 

Data  for-  Sampling  Plan  V:  N  =  [N^,  . ..,  N  ]  or  Y  =  [Y^,  . ..,  Yc] 

where,  Kr  =  [N^W^),  »r(W*k)],  and  Yr  =  [Yrl,  ....  YrR] 

Y  =  N(W*  )  -  N(W*  ,)  and  N  (0)  =  0  for  r  =  1,  c  and 

rs  rs  r,s-l  r  ’ 

s  =  1,  2,  k.  Here,  {W*p},.  •••»  {W*  }  are  waiting  times  of  independent 

NHPP's  each  with  mean  function  pQ(  • )  and  independent  of  N^t),  ...»  N  (t). 


Based  on  data  of  the  type  above,  one  finds  several  detection  rules 
to  be  related. 


Detection  Rules  D.I.l,  D.III.l  and  D.IV.l  Decide  N  +  S  iff 

s*  =  S(s"-  n"]2  [N]_1  >  X2  ( l-° ) . 

j  r  c~l 

For  Sampling  Plan  I,  III  and  IV,  n"  =  N(T*),  N(k a)  and  N(t,  ), 

r  K 

—  -i  C  ii 

respectively;  and  N  =  c  z  in  each  case. 

For  the  Sampling  Plan  II,  one  employs  Hartley's  (1950)  procedure 

p 

for  homoscedasticity,  which  is  applicable  here  since  2 )  ** 


Detection  Rule  D.II.l  Decide  N  +  S  iff 


[max  {vo(W  k)}]  [min  (u^Cw^)}]  >  h”,  the  appropriate  percentile  of 


o  rk 


Hartley's  distribution. 


The  final  detection  rule  is  related  to  determining  whether  or  not 
c  independent  negative  binomial  random  variables  have  distributions 
with  the  same  probability  p  (of  success). 

Based  on  the  likelihood-ratio  criterion,  one  derives 

Detection  Rule  D.V.l  Decide  N  +  S  iff  L*  >  x^Jl-®),  where 

—  .  c  z 

L*  =  2ck  ln(k+z)  -  2cz  In  ( — - — )  +  2  E  z  ln( - — — ) 

k  +  z  r=l  r  k  +  zr 

«  _  _i  c 

zr  =  Nr(wrk);  and  z  =  c  E 

One  other  detection  rule  is  of  some  interest  here. 

Meelis  (197U)  considers  a  situation  in  which  the  N  +  S  case  involves 
the  {6r }  constituting  a  random  sample  from  am  appropriate  continuous 
distribution.  Based  on  his  work,  one  arrives  at 

c  2  # 

Detection  Rule  D.V.2  Decide  N  +  S  iff  E  [N  (W*,  )]  >  d  , 

r=l 

c  *  *  # 

where  m  =  E  Nr(Wrk)  and  d  =  d  (a,  c,  m)  (the  appropriate  percentile 
of  Meelis'  distribution). 
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APPENDIX  -  Numerical  Examples 


[The  calculations  of  the  appendix  were  done  by  A.  Mason.] 

This  appendix  makes  use  of  seven  data  sets.  Those  designated  Data 
Set  1,  2,  3,  4  and  5  are  simulated  sets  of  data  from  NHPP's  with 
the  indicated  mean  functions.  The  last  two  data  sets  are  records  of 
seizures  for  two  different  epilepsy  patients. 

The  examples  are  grouped  by  problem  (designated  by  letters  A,  B, 
C,  D  and  E)  and  by  sampling  plan  (Roman  numerals  I,  II,  III,  IV  and 
V). 

For  example,  "Case  D.III."  refers  to  Detection  Problem  D  and 


Sampling  Plan  III. 


1 


*  ✓* 


Data 

Set  1 

Data  Set  2 

Data  Set  3 

Data  Set  4 

Data  Set  5 

i 

y(t) 

2 

=2. St 

y(t)=t2 

U(t)=4t 

y(t)=0.5t2 

U(t)=0.25t 

m 

.v. 

1 

0.6621 

1.5370 

0.2640 

1.6407 

1 . 8609 

2 

0.8156 

1.8415 

0.3535 

3.1687 

4.2283 

3 

1.1033 

2.2670 

0.8270 

3.1843 

4.9514 

4 

1.1425 

2.4329 

1.43] 3 

3.7095 

6.3765 

■ 

5 

1.3875 

2.6718 

1.7190 

3.7499 

7.4145 

6 

1.4569 

2.9192 

2.3334 

3.8316 

7.6608 

7 

1.6535 

2.9655 

2.7005 

4.4190 

8.8394 

8 

1.6554 

3.3051 

2.9981 

4.7912 

10.0257 

9 

1.7187 

3.3359 

3.0349 

4.8704 

10.7584 

10 

1.8768 

3.6561 

3.4578 

5.1042 

12.6299 

P 

11 

2.1242 

3.6623 

3.8193 

5.1537 

12.8330 

12 

2.1813 

3.8465 

4.1066 

5.3270 

12.9096 

\ 

13 

2.2706 

3.8716 

4.2841 

5.6595 

14.0877 

14 

2.3851 

3.9012 

4.2980 

5.8857 

14.9650 

15 

2.4161 

4.4260 

4.3277 

5.8861 

15.1432 

'  . 

16 

2.5016 

4.4830 

4.4119 

6.4762 

15.5248 

i 

17 

2.5161 

4.6552 

4.9939 

7.0334 

15.9582 

18 

2.5183 

4.7408 

5.7613 

7.0720 

16.2997 

19 

2.5231 

4.8441 

6.1637 

7.2839 

17.2403 

20 

2.5745 

4.8968 

6.4948 

7 . 3309 

17.3774 

21 

2.7927 

4.9904 

6.6561 

7.6767 

17.6973 

s> 

22 

2.9971 

5.0636 

7.1748 

7.7333 

18.6457 

■ 

23 

3.0902 

5.0779 

7.2400 

8.2574 

18.5428 

24 

3.1346 

5.0933 

7.3192 

8.5161 

18.5877 

25 

3.1652 

5.2479 

7.4199 

8.3744 

18.8340 

y; 

26 

3.1985 

5.4329 

7.7230 

9.0870 

18.8848 

■*.* 

27 

3.2867 

5.4493 

7.7686 

9.1491 

19.1638 

& 

28 

3.3588 

5.4559 

7,9864 

9.2653 

19.3535 

■ 

29 

3.4849 

5.5021 

8.3356 

9.7845 

19.7537 

30 

3.5187 

5.6157 

9.4819 

9.9437 

20.0240 

$ 

31 

3.5476 

5.6721 

9.7305 

10.0824 

21.1768 

32 

3.5925 

5.7419 

10.1198 

10.2001 

21.6865 

%■. 

V, 

33 

3.8160 

5.9278 

10.3810 

10.4780 

21.7489 

34 

3.8490 

5.9652 

10.5817 

10.6291 

22.2156 

■ 

35 

3.9115 

6.0453 

10.8142 

10.9366 

22.5020 

36 

3.9896 

6.1611 

10.9833 

11.4502 

23.1747 

■4'. 

37 

4.0754 

6.2578 

11.1056 

11.6886 

24.3404 

38 

4.1004 

6.3962 

11.3623 

11.8080 

23.7653 

39 

4.1185 

6.3982 

11.6609 

11.9752 

24.7081 

40 

4.1229 

6.4602 

11.7117 

12.0659 

25.1715 

j 

41 

4.2126 

6.5531 

12.1708 

12.2152 

25.2280 

42 

4.2949 

6.5856 

12.9102 

12.2362 

25.3121 

43 

4.2966 

6.5952 

12.9625 

12.3179 

25.3591 

»* . 
ji»  * 

44 

4.3525 

6.6012 

13.6038 

12.3855 

25.7556 

yi 

45 

4.4296 

6.6321 

14.5538 

12.4373 

25.7794 

46 

4.4329 

6.7129 

15.2888 

12.4553 

26.0114 

t*  47 

4.4472 

6.7171 

15.3994 

12.4829 

26.3572 

F 

'  48 

4.5433 

6.9320 

15.4766 

12.9268 

26,6164 

P: 

49 

4.5572 

7.0581 

15.5336 

13.0550 

26.9873 

>  50 

4.5681 

7.1849 

15.5897 

13.2268 

27.1501 

»*'  * 

Problem  A:  Class-fit,  or  Goodness -of- fit  with  Nuisance  Parameter. 


PN:  L  e  ^(yQ)  vs.  N  +  S:  L  t  ft(y0) 

2 

where  y^ft)  =  t  . 

Case  A. I :  Observe  N(t)  until  time  T*.  The  data  is  of  the  form 
Wj,  W2,  . ...  Wk  where  k  =  N(T*) . 


Detector  Statistics:  (i)  Kolmogorov- Smirnov 


D.  =  sup 
0<z<l 


i  k 

^  I  e  [z 

j-1 


y0(T*) 


■]  -  z 


^  K  -  S(k) 
PN 


Detection  Rules:  (i)  (A.I.l.)  Decide  N  +  S  iff  D^  >  \  a 

(ii)  (A. 1.2.)  Decide  N+S  iff  B  >  x?>  1  -  a 


Decide 


Decide 


Case  A. II :  Observe  N(t)  until  Nth  event.  The  data  is  of  the 


form  Wp  W2,  ....  WN« 


Detector  Statistics: 

(i)  Kolmogc^ov-Smirnov 


D 


N-l 


sup 

0<z<l 


1 

N-l 


N-l 

l 


j=l 


e  [z 


VWj> 

W 


z 


-v  K-S(N-l) 
PN 


(ii)  Neyman-Barton 

Bp  =  (N-l) [V  -  1/2,  V2  -  1/3, 


—  v-l 

VP  -  1/p+l]  l 


V  -  1/2 


V2  -  1/3 
VP  -  1/p+l 


P 


where 


N-l 

l 

j=l 


vv]r 

WJ 


and 


l  -  Cors) 


°rs  (r+s+1) (r+1) (s+1) 


Detection  Rules: 


(i) 

(A.II.l.) 

Decide 

N  .  S  iff  DfJ_I  >  dN.l  a 

(ii) 

(A. II. 2.) 

Decide 

N.S  iff  Bp  >  *^1(j 

For 

Data  Set  1 

(N  =  50)  the  following  values  are 

!  obtained: 

(i) 

D49  =  0.06 

Decide 

PN,  since  critical  value 

is 

d49,.01  =  0,23 

(ii) 

B2  *  0.60 

Decide 

PN,  since  critical  value 

is 

X2,.99  ’  °'21 

=  1 


2 

B,  =  1.16  Decide  PN,  since  critical  value  is  X_ 


nn 


z  ?<? 


For  Data  Set  3  (N  =  50)  the  following  values  are  obtained  (same 
critical  values  as  above) 


(i) 

°49 

=  0.04 

Decide 

PN. 

(ii) 

B2  = 

28.49 

Decide 

N  +  S. 

B,  = 

29.86 

Decide 

N  +  S. 

Case  A. Ill :  Regular  Sampling 

Observe  N(t)  at  regular  intervals  of  time  A,  2A,  kA 


Data  Set  1  (A  =  0.50) 

N(A)  =  0 

N(2A)  =  2 

N(3A)  =6 

N(4A)  =  10 

N(5A)  =  15 

N(6A)  =  22 

N(7A)  =  29 

N(8A)  =  36 

N(9A)  =  47 

k 

Detector  Statistic:  T  =  £ 

j=l 


Data  Set  2  (A  =  1.00) 

N(A)  =  0 

N(2A)  =  2 

N(3A)  =  7 

N(4A)  =  14 

N(5A)  =  21 

N(6A)  =  34 

N(7A)  =  48 


b.N(kA)^2 

yj  '  yQ(kA)  /  %  2 

b.N(kA)  PN  4-1 

U0(kA) 


where  y.  =  N(jA)  -  N([j  -  1]A)  j  =  1,  ...,  k 


and  b,  =  un(jA)  -  y_([i  -  1 ] A) 


Detection  Rule:  Decide  N+S  iff  T  >  X,  ,, 

-  k-l,l-a 

For  Data  Set  1  (k  =  9,  A  =  0.50)  we  compute  T  =  1.35.  Since  the 
2 

critical  value  is  X.  no  =  20.1  one  decides  PN. 

Of  •  y y 

Case  A. IV:  Equal -Distance  Sampling 

Observe  N(t)  at  times  t^,  t • ...  t^  where  the  t's  are  chosen 

so  that  u0(tj)  "  ^O^j  1^  are  eclual  for  every  j  =  1,  ...,  k;  tQ  = 

2 

If  we  assume  y^ft)  =  t  ,  a  convenient  choice  for  the  observation 

times  are  t^  =  for  j  =  1,  ...,  k. 

[y  (t.)  -  yn(t  ,)  =  1  for  all  j]. 


Data  Set  1  (k  =  10) 


Data  Set  2  (k  =  15) 


Data  Set  4  (k  =  15) 


Nftj)  =  2 

N(tp  =  9 

N(tj)  =  0 

N(t2)  =  5 

N(t2)  =  0 

N(t2)  =  0 

N(t3)  =  9 

N(t3)  =  1 

N(t3)  =  1 

N(t4)  =10 

N(t4)  =  2 

N(t4)  =  1 

N(ts)  =  12 

N(t5)  =  2 

N(ts)  =  1 

N(tfi)  =  IS 

N(t6)  =  4 

N(t6)  =  1 

N(t?)  =  20 

N(t?)  =  4 

N(ty)  =  i 

N(tg)  =  21 

N(tg)  =  5 

N(tg)  =  1 

N(tg)  =  22 

N{tg)  =  7 

N(tg)  =  1 

N(t10)  =  24 

N(ti0)  "  7 

N^t10^  =  1 

N(tn)  =  8 

N(tn)  =  3 

M(ti2)  =  9 

N(t12)  =  3 

N(ti3)  =  9 

N(tl3)  =  3 

N(tu)  =  11 

N(tu)  =  4 

N(tls)  =  12 

N(t15)  =  6 

Detector  Statistic:  T  = 

1_ 

Y 

l  (Y.  -  Y)2  ^  X2  . 
j=l  J  PN  Vl 

where 

Y.  =  N(t . )  -  N(t .  ,). 

J  1  3-1 

Detection  Rule:  Decide 

N 

+  s  iff  T  >  x*_1(  1.a 

For  Data  Set  1,  Y  =  2.4  and  we  get  T  =  6.77.  Since  the  critical 

2  — 
value  is  Xg  gg  =  21.67,  one  decides  PN.  For  Data  Set  2,  Y  =  0.80 

2 

and  we  get  T  =  10.45.  Since  the  critical  value  is  Xn  .  no  =  29.14, 


10 


one  decides  PN. 

Case  A. V. :  Same-Shape  Sampling. 

Observe  NHPP  N(t)  with  mean  function  SUgft)  at  t^mes  wf»  w2’  ’* 

w£,  where  the  w*’s  are  the  waiting  times  of  an  independent  NHPP 

N*(t)  with  mean  function  y^ft).  F°r  these  examples.  Data  Set  1  is 

2 

observed  at  the  waiting  times  of  Data  Set  2  (y(t)  ■  t  ).  Also,  two 

2 

more  independent  processes  with  y(t)  =  t  were  generated  to  provide 
times  of  observation  for  Data  Sets  4  and  5. 

Data  Set  1  Data  Set  4  Data  Set  5 


<  Dm 


11 


s  (n.  -  kp(c.)  ^  2 

Detector  Statistic:  T  =  J  - - - - —  ^  X_  , 

i=l  kp(ci)  PN  s" 

where  =  number  of  y^'s  in  (one  of  s  mutually  exclusive 

class  covering  the  non-negative  integers) . 

2 

Detection  Rule:  Decide  N  +  S  iff  T  X  .  , 

-  s-l,  1  — ot 

If  we  take  s  =  5,  then  for  i  =  1,  5 

Cj  =  (5j  +  i  -  1:  j  =  0,  1,  2,  ...  }.  Thus,  Cj  =  {0,  5,  10,  15,  ...  }, 

C2  *  (1,  6,  11,  16,  ...  }  etc.  For  Data  Set  1  we  have  k  =  16  and 

N(w*  ) 

N(wJ6)  =  47.  The  MLE  of  6  is  §  =  — ^  =  2.94.  We  also  have 

*  (1  +  0)”1  =  0.254  and  p^)  *  p(q)1'1  (1-qV1  for  i  =  1,  ...»  5. 

The  pCCj)  and  the  kp(C^)  (with  k  =  16)  are  computed  as: 


PCCj) 

=  0.330 

kp(Cj) 

=  5.280 

r~\ 

CM 

u 

w 

<D« 

=  0.246 

kp(C2) 

»  3.936 

P(C3) 

=  0.184 

kp(C3) 

=  2.944 

PCC4) 

*  0.137 

kp(C4) 

=  2.192 

P(C5) 

*  0.102 

kp(C5) 

=  1.632 

Tnf  n^  are: 


V  ^ 


-  12  - 


The  computed  statistic  value  is  T  =  6.915  and  the  critical  value  is 
X42  gg  =  13.277.  Therefore,  decide  PN. 


Problem  B:  Goodness-of-fit. 


PN:  y(-)  =  (8=B0)  vs‘  N  +  S:  t  W’^CB^Bg) 


Case  B.I. :  Observe  N(t)  until  time  T*. 


Detector  Statistic:  N(T*)  ^  PqCBqPqCT*)) 

PN 

Detection  Rule:  Decide  N  +  S  iff  N(T*)  £  a,  or  N(T*) 
where  a^  and  are  appropriate  percentiles  of  the  Poisson  distribution. 
For  Data  Set  1,  consider  the  following  detection  problem: 


PN:  6=1,  vs.  N  +  S:  6^1 


If  T*  =  3.00,  then  we  observe  N(3.00)  =  22.  Also,  y(T*)  =  B^y^fT*)  = 
2 

(1)(3.00)  =  9.00.  From  Poisson  tables,  the  critical  values  are 

aj  =  2  and  a^  =  18  (with  FAR  a  =  .0115).  Thus,  one  decides  N  +  S. 


Case  B. II . :  Observe  N(t)  until  Nth  event. 


Detector  Statistic:  T  =  2B„y„(W1T)  'v  X„lt  . 

-  00  N  PN  2N 

Detection  Rule:  Decide  N  +  S  iff  T  <  b^  or  T  >  where  b^ 

2 

and  b2  are  appropriate  percentiles  of  the  X2N-distribution. 

2 

For  Data  Set  2,  assume  Pg(t)  =  t  • 

We  observe  W^g  =  7.1849  and  consider  the  following  three  detection 


problems : 


(i) 

PN: 

8=2  vs . 

N  +  S: 

6  /  2 

(ii) 

PN: 

8  =  1  vs. 

N  +  S: 

8  t  1 

(iii) 

PN: 

3=0.5  vs. 

N  +  S: 

6  t  0.5 

In 

each 

case  the  critical  values 

are 

b 

=  x2 

1  100,. 005 

67.328 

and 

.  b 

=  x2  = 

2  100,. 995 

140.169. 

then 

compute 

(i) 

T  = 

2(2) (7. 1844) 2  = 

206.49 

and  decide  N  +  S. 

(ii) 

T  = 

2(1) (7 . 1849) 2  = 

103.25 

and  decide  PN. 

(iii) 

T  = 

2(.05) (7. 1849) 2 

=  51.62 

and  decide  N  +  S 

Case  B. III.  :  Regular  Sampling  (See  data  on  page  7) 

.2 


Detector  Statistic : 


1=1  i  o 


DetectioT  Rule:  Decide  N  ♦  S  iff  T  > 


Vl-a 


For  Data  Set  1,  we  will  consider: 

PN:  6  =  2.5  vs.  N  +  S:  6  r*  2.5 


3  «?«? 


j 

Vo 

<VVo>2 

frrW 

Vo 

1 

0 

0.62S 

0.391 

0.625 

2 

2 

1.875 

0.125 

3 

4 

3.125 

0.766 

4 

4 

4.375 

0.141 

5 

S 

5.625 

0.391 

0.069 

6 

7 

6.875 

0.016 

0.002 

7 

7 

8.125 

1.266 

0.156 

8 

7 

9.375 

5.641 

0.602 

9 

11 

10.625 

0.141 

0.013 

T  =  1.811 

2 

Since  the  critical  value  is  Xg  =  2.088  one  decides  PN. 

Case  B. IV. ;  Equal -Distance  Sampling  (see  data  on  page  9). 

Detector  Statistic:  N(tk)  ^  P0(60y0(tk)) 

Detection  Rule:  Decide  N  +  S  iff  N(tk)  £  a^  or  N(tk)  a ^ 
where  a^  and  a ^  are  appropriate  percentiles  of  the  Poisson  distribution 
Suppose  one  is  interested  in  testing 

PN:  B  =  1  vs.  N  ♦  S:  B  t  1 


for  Data  Sets  1  and  2.  Taking  k  =  10,  one  has  N1(tin)  =  24 


15 


•v' 


and  N2(tlQ)  =  7.  Since 
10.  From  Poisson  tables, 
(with  FAR  a  =  .0063), 
PN  for  Data  Set  2. 


tj  =  Jj,  it  follows  that  ^O^k^  =  Cl)  C v^XO) 2  = 
the  critical  values  are  a^  =  2  and  a2  =  20 
so  one  decides  N  +  S  for  Data  Set  1  and 


Case  B.V. :  Same-Shape  Sampling  (see  data  on  page  10) 

Detector  Statistic:  N(w*)  NB(k,p)  where  p  =  [1  +  B  ] 

PN 

Detection  Rule:  Decide  N  +  S  iff  N(w*)  <_  or  N(w*)  _>  b2> 
are  appropriate  percentiles  of  the  negative  Binomial  distribution. 
Consider  the  following  problem: 


V 


•y 


PN:  .  B  *  1  vs.  N  ■+  S:  B  ?  1 

For  Data  Set  1,  we  observe  N(w*^)  =  47.  Thus  k  =  16  and  p  *  ^  *  -g-  =  y  . 

To  compute  critical  values,  we  must  find  b^  and  b2  such  that 

P(x  <_  bj)  =  y  and  P(x  b2)  =  y  where  X  *v>  NB(16,  y) .  A  corrected 

(Gram- Chari ier)  Poisson  approximation  will  be  used  (see  Johnson  and  Kotz, 

-v  bi  -  NP 

Discrete  Distributions,  pg.  129).  P(X  <_  b^)  ^  P(Y  £  b^) - 2^  V~P)  P^Y  =  bl^ 

where  N  =  16,  P  =  y-E —  =  1  and  Y  ^  Pq(NP).  Thus, 

.  b  -  16 

P(X  <  bx)  ^  P(Y  <  bx)  -  — -  P(Y  =  bx) 

where  Y  ^  PQ(16). 


From  Poisson  tables, 
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P(Y  =  b2  -  1) 

where  N,  P  and  Y  are  defined  as  above.  Thus, 

P(X<  28)  ^  .9978  -  3(.0019)  =  .9921 

so  b2  =  29  with  approximate  tail  probability  .0079. 

The  critical  values  are  therefore  b^  =  6  and  b2  =  29  with 
FAR  a  =  .0184.  Since  we  observed  N(w*6)  =  47,  decide  N  +  S. 

Problem  C:  Two-sample 

PN:  L  =  L-  vs.  N  +  S:  L.  /  L- 

1  C61=62)  Z(B^B2) 

Case  C . I . :  Observe  N(t)  until  time  T*. 

Detection  Statistic:  Conditionally  given  N(T*)  +  M(T*)  =  N, 

N(T*)  ^  B(N,  1/2)  (under  PN) 

Detection  Rule:  Decide  N  +  S  iff  N(T*)  <  c^  or  N(T*)  >  c2> 


P(X  <  6)  ^  .0040  -  ■  (.0026)  =  .0105 

so  b^  =  6  with  approximate  tail  probability  .0105. 

To  find  b2,  we  must  satisfy: 

P(X  >  b2)  =  1  -  P(X  <  b2  -  1)  =  ~ 

or  P(X  <  b2  -  1)  «  1  -  | 

(b2  -.1)  - 

P(X.  <  b2  -  1)  *  P(Y  <  b2  -  1) - 


can  be  determined  using  the  Normal  approximation  to  the  Binomial. 


cx  =  N(i/2)  -  za/2  rtUTzTUTz) 
c2  =  N  C 1/2)  +  z a/2  /N(l/2)(l/2) 

If  Data  Sets  1  and  2  are  observed  until  T*.=  3.00,  we  arrive  at 
N(T*)  =  22  and  M(T*)  =7.  For  a  =  .01 

Cj  =  29(1/2)  -  2.S76  /29(l/2) (1/2)  =  7.56 

c2  =  29(1/2)  +  2.576  /29(l/2) (1/2)  =  21.44 

Since  N(T*)  =  22,  one  decides  N  +  S. 

Case  C. II. :  Observe  N(t)  until  Nth  event. 

mW 

Detector  Statistic:  T  =  ^  F(2n,  2m) 

-  ny0(V  PN 


Detection  Rule:  Decide  N  +  S  iff  T  <  F^2n  2m 

T  >  F(2n,  2m,  l-a/2) ’ 

For  Data  Sets  1  and  2  assume 


or 


Ly  i2  E  n(yQ)  and  yQ(t)  =  t2. 


Taking  n  *  40  and  m  *  50  we  find: 


80,100, .995 


Since  the  critical  values  are  Fg(j  1Q0  0Q5  =  0.70  and  F 

decide  N  +  S. 

Case  C. Ill . :  Regular  Sampling 
Detector  Statistic:  Conditionally  given 

N(nAx)  +  M(mA2)  =  N,  NfnA^  n,  B(N,p)  (under  PN) 

V"V 

where  p  -  7  y0(mA2)  • 

Detection  Rule;  Decide  N  +  S  iff  NfnAp  <  Cj  or  NfnA^^)  >  c2 
where  c ^  and  c2  can  be  determined  using  the  normal  approximation 
to  the  Binomial: 

C1  =  Np  -  za/2  ^NpCl  -  p) 

c2  =  NP  +  za/2  ^Pt1  -  P) 

If  Data  Sets  1  and  2  are  observed  (with  n  =9,  A^  =  0.50,  m  = 
A2  =  1.00)  we  arrive  at  N(nA^)  =47  and  M(mA2)  =  48.  Thus, 

N  =  N(nAj)  +  M(mA2)  =  47  +  48  =  95  and 


A  C 


bj  =  (36) ( .4)  -  2.576  /36(.4)(.6)  =  6.83 

b2  =  (36) (.4)  +  2.576  /35(.4)(.6)  =  21.97 

Since  N(t^)  =  24,  one  decides  N  +  S. 

Case  C.V. :  Same-Shape  Sampling 

Detector  Statistic:  Conditionally  given 

N(W*)  +  M(V*)  =  k,  N(W*)  *  F* 
n  "  n  PN 

where 

PP*  {N(W*)  =-  s|N(W*)  +  M(V*)  =  k} 


Detection  Rule;  Decide  N  +  S  iff  N(W*)  <  c*  or  N(W*)  >  c*  wher 
c*  and  c*  are  the  largest  and  smallest  non-negative  integers  such  that 

c*-l 

ll  PF*[N(W*)  =  s|N(W*)  +  M(V*)  =  k]  <  | 

3=0 

k 

l  Pp*[N(W*)  =  slN(w£)  +  M(V*)  =  k]  if 

s=c*+l 

For  Data  Sets  4  and  5  we  observe  N(W*^)  =  5  and  M(V*^)  =2,  so 
k  =  N(W*J  +  M(V*£)  =7.  To  compute  critical  values  c*  and  c*: 


170544 

12620256 


.0135 


Thus,  for  FAR  a  =  .027,  the  critical  values  are  c*  =  1  and  c 
Since  N(WJ^)  =  5,  decide  PN. 

Problem  D:  Constant  Intensity  Rate. 

PN:  Bj  =  YqB2  vs.  N  +  S:  Bj  YQB2 

Case  D. I. :  Observe  N(t)  until  time  T*. 

Detector  Statistic:  Conditionally  given 

N(T*)  +  M(T*)  =  N,  N(T*)  ^  B(N,p)  (under  PN) 

Yo 

where  p  =  .  v 


Detection  Rule: 


where  c^  and  c2 


with  p  =  ■ 


yr 


Decide  N  +  S  iff  N(T*)  <  Cj  or  M(T*)  >  c2> 
are  determined  as  in  Case  C.I,  replacing  p  =  1/2 

2  5 

Suppose  we  take  Yq  =  2.5,  then  p  =  ~f "5 


Using  Data  Sets  1  and  2,  consider  the  following  problem: 


PN:  ux(t)  =  (2.5)y2(t)  vs.  N  +  S:  y^t)  /  (2.5)y2(t) 


As  before,  we  observe  N(T*)  =  22  and  M(T*)  =  7  for  T*  =  3.00. 


K>* 


The  critical  values  are  computed  below:  (a  =  .01). 


cl  =  29(0.71)  -  2.576  /29(0. 71) (0.29)  =  14.45 
c2  =  29(0.71)  +  2.576  /29 (0.71)  (.0.29)  =  26.98 

Since  N(T*)  =  22,  one  decides  PN. 


Case  D. II. :  Observe  N(t)  until  Nth  event. 


Detector  Statistic: 


"VW 

nU0CW*) 


*  F(2n, 
PN 


2m) 


As  in  Case  C.II,  we  observe  W4Q  =  4.1229  and  W*Q  =  7.1849.  If 
Yq  =  2.5,  the  computed  value  of  the  detector  statistic  is 


T  =  50(2.5)  (4.1229) 
40(7. 1849)2 


1.029 


Since  the  critical  values  are  Fg0>100> .005  =  0.70  and  F80j100).995 
decide  PN. 


Case  D.III. :  Regular  Sampling 
Detector  Statistic:  Conditionally  given 

N(nAp  +  M(mA2)  =  N,  N(nA^)  %  B(N,p)  (under  PN) 


where 


YoVnV 

Y0^0^n^l^  + 


Detection  Rule:  Decide  N  +  S  iff  N(nAp  <  ^  or  N(nAp  >  c. 
where  Cj  and  c^  are  computed  as  in  Case  C.III. 

For  Data  Sets  1  and  2  we  consider: 


PN:  =  (2.5)62  vs.  N  +  S:  ^  (2.5)82 

2 

Again,  N  =  95  and  p  =  - -  =  0.51. 

(2.5) (4.5 j  +  (7.0) ^ 

Thus,  (for  a  =  .01) 

Cj  *  (95) (.51)  -  (2.576)  /95(.51) (.49)  *  35.75 
c2  =  (95) ( .51)  +  (2.576)  /95( .51) ( .49)  =  60.80 
Since  N(nAp  =  47,  one  decides  PN. 


Case  D. IV. :  Equal -Distance  Sampling 


Detector  Statistics:  Conditionally  given 

N(t  )  +  M(t  )  =  N,  N(t  )  'v  B(N,p) 
n  m  n 


WV 


Wo<V  *  W 


where 


24 


Detection  Rule:  Decide  N  +  S  iff  N(t  )  <  b,  or  N(t  )  >  b- 
-  n  1  n  2 

where  and  b2  are  determined  as  in  Case  C.IV.  Suppose  one  wants 

to  test 

PN:  Bj  =  2.502  vs.  N  +  S:  6^  2.502 

for  Data  Sets  1  and  2.  For  n  =  10  and  m  =  15,  we  observe 

N(tio)  =  24  and  M(t^,-)  =  12.  Since  Ug(t)  =  1  ant^  tj  =  ^ '  we  ^ave 


(2.S)(^0)2+(v^5)2 


(2.5) (10) 
(2.5) (10)+15 


=  0.63 


Thus  (for  a  -  .01) , 

bl  =  36( .63)  -  2.576  /36( .63) ( . 38)  =  15.02 
b2  =  36 (.63)  +  2.S76  /36( .63) (. 38)  =  29.98 
Since  N(t1Q)  *  24,  one  decides  PN. 

Case  D.V. :  Same -Shape  Sampling 

Here,  two  methods  are  available:  Method  1,  based  on  the  likelihood 
ratio  and  Method  2,  using  a  normal  approximation  to  the  Binomial  distri¬ 
bution.  Sampling  Plan  V  is  slightly  modified  for  Method  2.  Instead  of 
generating  two  independent  ’’noise"  processes  to  provide  times  of  observa¬ 
tions,  we  observe  the  two  processes  of  interest  at  the  waiting  time  of  a 
single  independent  same-shape  process. 


Consider  the 


27 


N(w*)  ^  B(N,  y— - 
K  PN  1  Y< 


-)  (for  large  N) 


Detection  Rule:  Decide  N  +  S  iff  Z  >  z  or  Z  <  -z 

-  a/ 2  a/2 


Again,  consider  the  problem. 


PN:  8.  =  58.  vs.  N  +  S:  8.  t  56. 

14  14 


This  time  we  observe  both  Data  Sets  1  and  4  at  the  waiting  times 


from  Data  Set  2.  If  we  take  k  =  16,  wjg  =  4.4830  and  our  final 


observations  are 


N(wJ6)  =  47  and  M(w*6)  =  7. 


Thus, 


47  -  5(7) 


=  0.731 


'(47  +  7)5 


and  we  decide  PN. 


Consider  two  epileptics  (patients  1  and  2)  who  have  seizure  times 


given  on  pages  2-3.  Assume  the  processes  can  be  modeled  as  NHPP’s 


with  mean  function  of  the  form 


^(t)  =  B.tv 


i  =  1,2 


We  want  to  test  the  following  situation: 


PN:  6X  =  0.2582  vs.  N  +  S:  Bj  *  0.2562 


An  independent  NHPP  with  mean  function  y(t)  =  t0'5  is  generated  to 
provide  times  of  observation.  For  k  =  50,  we  have 

w*q  =  239.40  min.  =  4.0  hours 

and  we  observe 


Problem  E:  c-sample 


PN:  Lj  *  ...  ■  vs.  N  +  S:  L^'s  not  all  equal 

Case  E. I. :  Observe  N(t)  until  time  T*. 

1 

where  N  *-  T  N.(T*). 
c  j=l  > 

2 

Detection  Rule:  Decide  N  +  S  iff  T  >  X  ,  , 

-  c-1, 1-a 

Consider  the  case  c  =  4,  where  the  four  samples  consist  of  the 
points  of  Data  Set  3  observed  at  time  0  -  4.0,  4.0  -  8.0,  8.0  -  12. 


Detector  Statistic: 


c 

T-  l 

j=l 


(N. (T*)  -  N)2 


12.0  -  16.0.  Thus  T*  =  4.0  and  we  observe  N^T*)  =  11,  N2(T*)  =  17, 

N3(T*)  =  12  and  N4(T*)  =  10.  (N  =  12.5).  The  computed  statistic 

2 

value  is  T  =  2.32.  Since  the  critical  value  is  X,  QQ  =  11.35,  one 
decides  PN. 


Case  E. II. :  Observe  N(t)  until  Nth  event. 
Detector  Statistic:  T  = 


max  {yQ(WiN):i  =  1,.. 


min  {uQ(WiN)':i  =  1,.. 


.  c} 
.  c> 


Detection  Rule 
100(1  -  a)%  point  of  the 


Decide  N  +  S  iff  T  >  b*  where  b*  is  the 
S2  max 


2 

S  min 


distribution  (Hartley) . 


If  Data  Set  3  is  divided  into  4  samples  of  12  observations  each,  we 
have 


Wl,  12 
w2, 12 
W3,12 


4.1066 

3.2126 

3.6641 

4.4933 


Since  L ^  e  Q(yQ) ,  i  =  1,  ...,  4  with  yQ(t)  =  t  we  compute 


T  = 


4.4933 

3.2126 


1 . 3986 


From  tables  of  the 


S2  max 
S2  min 


distribution,  a  critical  value  of 


7.6  is  obtained  (for  a  =  .01'  and  11  degrees  of  freedom). 


-  30  - 


We  therefore  decide  PN. 


Case  E. III. :  Regular  Sampling 


Detector  Statistic:  T*  =  £  T.  ^  X2  (ck  -  k) 

j=l  3  PN 


c  (N.(jA)-N.)2  _  ,  c 

where  T.  *  I  - — — — * -  and  N.  =  —  £  N.(jA) 

3  i=l  N.  3  c  i=l  1 


Detection  Rule:  Decide  N  +  S  iff  T*  >  X*k  ^  j  a 

If,  for  Data  Set  3,  we  take  A  =  0.50  and  k  =  10,  N^jA),  N2(jA) 
and  N3(jA)  can  be  generated  conSectutively: 


Nj  ( j  A) 


N2(jA) 


N3(jA) 


0.50 

2 

0 

2 

1.33 

2.00 

1.00 

3 

1 

5 

3.00 

2.67 

1.50 

4 

3 

7 

4.67 

1.86 

2.00 

5 

4 

9 

6.00 

2.33 

2.50 

6 

8 

10 

8.00 

1.00 

3.00 

8 

11 

12 

10.33 

0.84 

3.50 

10 

12 

12 

11.33 

0.24 

4.00 

11 

12 

13 

12.00 

0.17 

4.50 

16 

13 

13 

14.00 

0.43 

5.00 

17 

14 

14 

15.00 

0.40 

T*  =  11.94 

2 

Since  the  critical  value  is  X  2Q  gg  =  37. 57  one  decides  PN. 
Case  E.IV. :  Equal -Distance  Sampling 
Detector  Statistic:  W  =  -2  log  1 


l  VV 

W  l0* 


c  N. (tv) 

■]  *2  l  N.(tk)  logt-i^S-] 


.  2 

Detection  Rule:  Decide  N  +  S  iff  W  >  X  .  , 
-  c- 1 , 1-a 

Suppose  we  take  k  -  10  and  c  =  3  with  Data  Sets  1,  2  and  4 
consideration.  Then  we  have 

NjCt^)  =  24,  N2(t1Q)  =  7  and  Nj(t1(j)  =  1 

The  statistic  is  computed  as  follows: 

W  =  -2  [32  log(|§0  ]  +  2  [24  log(fi)  +  7  logf^) 

+  1  logfji-)]  =  28.29 
2 

Since  the  critical  value  is  X2  gg  =  9.21,  one  decides  N  +  S. 
Case  E.V. :  Same-Shape  Sampling 

f\j  7 

Detector  Statistic:  L  =  -2  log  1  ^  X  . 

PN  C-1 


under 


where 


2  <?<? 


_  Y  e  X. 

L  =  2ck  log(k  +  X)  -  2cX  logC-^  +2  [  X,  log (-—!—) 

k+X  i=l  1  k  i 


c 


-  2k  l  log (k  +  X.) 
i=l  1 


and  X.  =  N.(w*). 
l  l  k 

2 

Detection  Rule:  Decide  N  +  S  iff  L  >  X  ,  , 

- .  . . . .  ■  c- 1 , 1  -a 

Suppose  we  consider  Data  Sets  1,  4,  and  5  with  k  =  16.  We  then 

have 

X1  3  Nl(w165  3  47  X  =  18.0 

X2  -  N2(w*6)  =5 

X  ■=  N,(w*  )  »  2  cX  log( — - — )  =  541og(0.53)  =  -34.34 

*  *  AD  k+X 


L  =  2(3) (16)  log (34.0)  -  2(-34.34) 

+  2 [47  log  (||)  +  5  log  (^f)  +  2  log  (y|)]  -  2(16) [log  63  +  log  21  +  log  18] 

=  338.53  +  68.69  -  50.68  -  322.50  =  34.04. 

2 

Since  the  critical  value  is  X2  gs  =  5.991  decide  N  +  S. 


